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Abstract 

In this work, we consider the random access of nodes adapting their transmission probabihty based on the local 
channel state information (CSI) in a decentralized manner, which is called CARA. The CSI is not directly available 
to each node but estimated with some errors in our scenario. Thus, the impact of imperfect CSI on the performance 
of CARA is our main concern. Specifically, an exact stability analysis is carried out when a pair of bursty sources are 
competing for a common receiver and, thereby, have interdependent services. The analysis also takes into account 
the compound effects of the multipacket reception (MPR) capability at the receiver The contributions in this paper 
are twofold: first, we obtain the exact stability region of CARA in the presence of channel estimation errors; such 
an assessment is necessary as the errors in channel estimation are inevitable in the practical situation. Secondly, 
we compare the performance of CARA to that achieved by the class of stationary scheduling policies that make 
decisions in a centraUzed manner based on the CSI feedback. It is shown that the stability region of CARA is not 
necessarily a subset of that of centralized schedulers as the MPR capability improves. 

Index Terms 

Channel-aware random access, channel estimation errors, stability region, multipacket reception capacity 

I. Introduction 

INCREASING demand for high data rate to support a wide range of services in wireless data networks has led 
to the exploitation of diversity amongst users lH]. The diversity gain arises from the fact that wireless links 
experience random fading due to the constructive/destructive effect of multipath signal propagation and, thereby, 
there is always a user having better channel quality than the others at any time ||2l. A downUnk scheduler exploiting 
such diversity gain is called the opportunistic scheduler O-IH. Similar concept can be applied to the uplink 
communication but it is needed to have a centralized controller gathering channel state infomiation (CSI) from 
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distributed users, making a centralized decision, and distributing the decision information back to the distributed 
users. 

The ALOHA protocol, the simple scheme of attempting transmission randomly, independently, distributively, and 
based on simple ACK/NACK feedback from the receiver, has gained continued popularity in distributed multiaccess 
communication systems since its creation [6|. It is mainly due to its simplicity and the fact that it does not require 
centralized controllers. It also provided a basis for following enhanced schemes such as the Carrier Sense Multiple 
Access (CSMA) that combined ideas of fixed allocation with reservations via contention through random access. 
However, the random access systems have been built upon simplistic assumptions on the physical-layer such as the 
collision channel model, and the diversity in channel quality across different users was mostly out of consideration. 
Noticeably, there is a recent line of work on exploiting CSI under random access framework, which is called 
channel-aware random access (CARA) fTl- lfTOl . The CARA allows the distributed nodes to adjust their random 
access probability based on the local CSI. 

To put our contribution in perspective, we start with some background on CARA. In |7|, it is assumed that 
each user has perfect local CSI and transmits only when the channel gain exceeds a certain threshold. The 
main contribution of the work is the characterization of the throughput scaling law for the system with infinitely 
backlogged users, i.e., users have packets to transmit at any time. Limitations of this work include the collision 
assumption made for the analysis when two or more users transmit at the same time. This assumption cannot be 
validated in general wireless communications since a transmission may succeed even in the presence of interference 
|[m - |[T4l . Also users were assumed to be symmetric in channel statistics. In [Si, a similar problem with that in Q 
was considered but with additional multipacket reception (MPR) capability ifTSl - lfTTl . The advent of the multiuser 
detection technique for separating signals from the superposition of multiple received signals enables a receiver 
to correctly decode more than one packets simultaneously transmitted from different users. However, like [7 |, this 
work is based on the assumptions that users are always backlogged and symmetric in channel statistics. Note that in 
the system with bursty input traffic, it is not straightforward to know how users interfere with each other since they 
transmit only when having non-empty queues, which is more reahstic than assuming always backlogged queues. 
In 191, the stabihty regioiu of the system comprised of users having bursty packet arrivals and asymmetric channel 
statistics was obtained under the collision assumption and further extended in |[TOl to the case with MPR capability. 
However, the analysis in both ||9l and ifTOl is limited to the two-user scenario. This limitation is mainly due to the 
complex interaction between network queues. More specifically, the service process of individual queues depends 
on the status of the others. This is why most previous work on the stability of interacting systems has been focused 
on small-sized network and only bounds or approximations are known for networks with an arbitrary number of 
users lfT8l-ll24l. 

In this work, we focus attention on the impact of imperfect CSI on the stability region of CARA for bursty 



'stability region is a set of arrival rates maintaining the network queues in a bounded region. A formal definition is given in Section HH 
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input traffic. It needs to be empliasized that tiie analysis in all the above-mentioned work on CARA was performed 
based on the ideal assumption that the perfect local CSI is available at each user fTl- lfTOl . In reality, however, 
the CSI is obtained through an estimation and any kind of estimation is imperfect as long as there is randomness 
in the observed signal. Consequently, the occurrence of errors in estimation is inevitable and the performance of 
CARA would certainly depend on the accuracy of channel estimation. To see the effect of imperfect CSI, we allow 
channel estimation errors in the two-state time-varying channel model considered in ||9], flO \. in which each user i 
transmits with probability pi when its channel state is good and the packet queue is non-empty. Note that the errors 
either deprives the chance to utilize the good channel state when falsely estimated to be bad or causes unnecessary 
interference to the other user when its actual channel state is bad but falsely estimated to be good. There is some 
related work on CARA with imperfect CSI; in ||25l . the two-user system was considered under the collision channel 
model and, in |26], a system with an arbitrary number of users was considered but with always backlogged queues 
and symmetric channel statistics. 

Our contribution in this work can be summarized as follows. We first introduce the realistic effect of practical 
channel estimation into the stability analysis of CARA. The analysis also takes into account the compound effect 
of the MPR capability, which depends not only on the set of transmitters but also on their instantaneous channel 
states. The derived stability region describes the theoretical limit on rates that can be pushed into the system 
while maintaining the queues stable at given channel estimation error rates and MPR probabilities. Secondly, by 
comparing with the case of having perfect CSI, we identify the loss due to the imperfect CSI on the stability region 
of CARA. Finally, the stability region of the longest connected queue (LCQ) policy ITTl . which is a throughput- 
optimal policy that can stabilize the system whenever the stability is attainable, is derived again in the presence 
of channel estimation errors. The LCQ policy schedules a user having longest queue among those whose channel 
is connected and, thus, requires queue length and channel state information feedback to the centralized controller. 
Interestingly, we observed that the stability region of CARA, a fully distributed policy, is not always a proper subset 
of that of the LCQ policy. This is when relatively strong MPR capability presents. 

The rest of the paper is organized as follows. In Section |II] we present the system model and revisit the notion 
of stability. In Section |llll we describe our main result on the stability region of CARA in the presence of channel 
estimation errors. The proof of our main result is presented in Section |IV] which is based on the stochastic dominance 
technique previously introduced in |[T9l to deal with interacting queues. In Section |Vl the stability region of CARA 
is compared to that of the LCQ policy. Finally, we draw some conclusions in Section |Vll 

II. System Model 

We consider a multi-access system consisting of N nodes and a common receiver. The stability analysis of 
CARA is done for = 2 as illustrated in Fig. [T] whereas N is set to an arbitrary positive integer for the analysis 
of the LCQ policy in Section |Vl Each node i has an infinite size queue for storing the arriving packets that have 
fixed length. Time is slotted and the slot duration is equal to one packet transmission time. Let Qi{n) denote the 
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number of packets buffered at Si at the beginning of the n-th slot which evolves according to 

Qi{n + 1) = max[Qi{n) - fii{n), 0] + Ai{n) (1) 

where the stochastic processes {Ai^(?^)}^o ^^'^ {^ii^)}'?^=o sequences of random variables representing the 
number of arrivals and services at Sj during time slot n, respectively. The arrival process {Ai{n)}'^^Q is modeled 
as an independent and identically distributed (i.i.d.) Bernoulli process with E[Ai{n)] = Aj, and the processes at 
different nodes are assumed to be independent of each other. The service process {Aij(n.)}J^Q depends jointly on 
the transmission protocol and the underlying channel model, which governs the success of transmissions. 

The channel between node i and the receiver is randomly time-varying and its state at time slot n is denoted by 
Cj(n) and let C(n) = {Cj(n), . . . , CN{n)}. We assume that channels hold their state for the duration of a slot and 
potentially change on the slot boundarie^ As in the previous work |j9], lITOl . we model the time-varying channel as 
a discrete-time stochastic process taking values from {G, B}, which represents good and bad states, respectively. 
The channel processes at different nodes are assumed to be independent of each other but the realizations of a 
channel at a node at different time slots are not necessarily independent. If fact, a channel process can be arbitrarily 
correlated over time as long as stationary. We denote by vrp and 'Kf{= 1 — vrp) the steady-state probabilities that 
channel i is in good and bad states, respectively. The transmission control policy considered in [10] is studied 
again in which node i transmits with probability when Ci{n) = G and its queue is non-empty. We denoted by 
Ci{n) the estimated channel state over the link between node i and the receiver, and node i is now transmitting 
with probability pi if Ci{n) = G and its queue is non-empty. It is obvious that the performance of this adaptation 
would highly depend on the accuracy of the channel estimation. Let us define ef = Pr[C'i(n) = B\Ci{n) = G] 



^This assumption can be validated when the channel coherence time is relatively longer than the slot duration. 
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and ef = Pr[Cj(n) = G\Ci{n) = B], which are the probabiUties of falsely estimating the channel state and let 
ef = 1 - ef and ef = 1 - ef . 

The success of a transmission depends on the underlying channel model. The MPR model used in this work 
enables the probabilistic reception of simultaneously transmitted packets owing to the multi-user detection ||28l . 
Unlike the static MPR model used in ifTSl - lfTTl . in which the time- varying nature of wireless channels was 
disregarded, our model captures not only the effect of interference but also the instantaneous channel states of 
the transmitters. Denote by A/tx and Ctxi'n) the set of transmitting nodes and their channel states, respectively. 
Then, the probability that a packet transmitted from node i is successfully decoded at the destination is described 
by 

where 7i|Ctx(n) is the signal-to-interference-plus-noise-ratio (SINR) of the signal transmitted from node i at the 
receiver given the channel states of the transmitters Ctxin) and the threshold for the successful decoding 6, which 
depends on the modulation scheme, target bit-error-rate, and the number of bits in the packet |29|. We assume 
throughout the paper that the success probability when a node's own channel state is bad is negligible due to such 
as the deep fading regardless of the interference. This is the situation when the signal-to-noise-ratio (SNR) itself is 
below the threshold 6 for the successful decoding. Thus, for = 2, we are particularly interested in the following 
reception probabilitiej^ 

Ql\{G}, Ql\{G,B}, Ql\{G,G}, Q2\{G}^ Q2\{B,G}^ Q2\{G,G} 

Note that the success probability when a signal is transmitted in the presence of interference cannot exceed the 
probability when it is transmitted alone. Moreover, since we are considering the multi-access to a common receiver, 
the fact that a node's channel is in the good state implies potentially higher interference level to the other node 
if they transmit at the same time. Therefore, the following relations hold: qi\{G} > Qi\{g,b} > Qi\{G,G} ^^id 

92|{G} > Q2\{B,G} > Q2\{G,G}- 

We adopt the notion of stability used in EH in which the stability of a queue is equivalent to the existence of 
a proper limiting distribution. In other words, a queue is said to be stable if 

lim Pr[Qj(n) < x] = F{x) and lim F{x) = 1 

n— >oo x~>oo 

Loynes' theorem, as it relates to stability, plays a central role in our approach |]30]. It states that if the arrival 
and service processes of a queue are strictly jointly stationary and the average arrival rate is less than the average 
service rate, the queue is stable. For given channel estimation error rates e and transmission probabilities p of the 
nodes, the stability region &{e,p) is defined as the set of arrival rate vectors A = (Ai, A2) for which all queues in 
the system are stable. The stability region of the system G(e) is defined as the closure of &{e,p) over all possible 
transmission probability vectors, i.e., ©(e) = Upe[o i]^ 



^For example, q2\{B,G} is the probability that the transmission by node 2 is successful when Ci(n) = B and C2{n) = G. 
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III. CARA WITH Imperfect CSI 

This section describes the stability region of CARA in the presence of channel estimation errors. As noted earlier, 
the service process of a queue depends on the status of the other, which makes the analysis challenging. The proof 
of the results, which is based on the stochastic dominance technique [19], is presented in the next section. Let us 
define 

*i = ^f'ef {qi\{G} - Ql\{G,G}) + - Ql\{G,B}) 

% = {Q2\{G} - <12\{G,G}) + Trf ef {q2\{G} - <12\{B,G}) 

which are shorthand notations to simplify the description of our main results. 

Lemma 3.1: The stability region &{e,p) of CARA at given channel estimation error rate vector e and transmis- 
sion probability vector p is the union of the following subregions: 



TZi = I (Ai, A2) : Ai < vrfefpi Iq^^G} " ^fif(^^^^^^3^|^ I ' < ^2 4p2 {q2\{G} " 



and 



n2 = |(Ai,A2) : Ai < vrfef pi {q^^^c} " ^1^2) , A2 < nfe§P2 ^g2|{G} " vrfef {q^i{G} - nP2) 

Proof: The proof is presented in the next section. 
Let us define the following points in the two-dimensional Euclidean space: 



^'^ Q2\{G}iQl\{G}-n? T^2^nq2\{G} \ 



Pi 



*29l|{G} 9l|{G} 



Trf ef%q^{G} ^2^e^9i|{G}fe|{G} " ^DM 
V 92|{G} *l92|{G} / 



P3 = (vrf fe|{G} - vr^ef feliG} " ^D) (4) 

which are aU in the first quadrant. 

Theorem 3.1: If -^^^ + ^^^^ — 1' boundary of the stability region ©(e) of CARA at a given channel 
estimation error rate vector e is described by three segments: (i) the line connecting Py = (0, 7r2^e2^g2|{G}) ^''^^ ^i' 
(ii) the curve 

\ ZG^G ^1 + \ rC^G ^2 = ,Jqi\{G}q2\{G} (5) 
V ^1 V 2 ^2 

from Pi to P2, and (iii) the line connecting P2 and Px = (Trf^e^giHG}, 0). If -^^^ + < !> it described 

by two lines: (i) the line connecting Py and P3 and (ii) the line connecting P3 and Px- 

Proof: The proof is presented in the next section. ■ 
Corollary 3.1: If — ^ — | — > 1, the stability region (3(e) is non-convex. If — ^ — I — < 1, it is a convex 

SlKG} <?2|{G} J b \ J <ll\{G} 92|{G} — 

polygon. Specifically, when ^^^-^ + = region becomes a right triangle. 
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This corollary can be easily verified by comparing the slopes of the lines from Py to Pi and from P2 to Px and 
those from Py to i-3 and from P3 to Px- Specifically, if ^^-^j^ + q^^^^iy — curve dSJl shrinks to a point whose 

coordinate is identically described by both Pi and P2 and the slopes of the lines from Py to Pi and from P2 to 
Px become identical. 

Consider the case when perfect CSI is available. This can be viewed as a special case of our model with 
e = 0, where is the vector of zeros. By substituting e = into Theorem 13.11 we can obtain the stability 
region for the case with perfect CSI, which reconfirms the previous results obtained in ||9] and HIOJ . For the 
comparison's sake, let us consider the case when CSI is not available and, hence, each node has to make decisions 
on transmission independent of the underlying channel states. This corresponds to the original ALOHA in which 
each node transmits with probability pi regardless of the underlying channel states whenever its queue is non-empty. 
Thus, at a given set of transmitters, the success probability of each node is given as a constant, which is obtained 
by taking the average over the stationary distribution of the channel states. Denote by qf and qf^ {i G {1,2}) the 
transmission success probabilities seen by node i when it transmits alone or along with the other node j(/ i). For 
the two-node case, it is obtained as ql = vrf gi||G}, qT = ■^iT^2<1i\{G,G} + ■^i'^iQi\{G,B}, qI = 7r^92|{G}> and 
= T^iT^2 12\{G,G} + '^i'^2^2\{B,G}- Also define Aj = — g™, which is assumed to be strictly positive without 
loss of generality. The following theorem obtained in |[T6l describes the stability region of the original ALOHA 
for the case with static MPR channels, but it is also applicable to the system with time-varying channels but when 
the CSI is unavailable. This is because the success probabilities are given as constants over time. We especially 
denote the stability region for this case by ©(0) since the notion of channel estimation errors for the case with no 
CSI is not valid. 

Theorem 3.2: If ^ + ^ > 1, the boundary of the stability region of ALOHA with no CSI, denoted by 6(0), 
is described by three segments: (i) the line connecting Py = (0, g|) and Pi = , ^^)> (ii) the curve 

^/A^+^/A^ = from Pi to P2 = ^^%^|^), and (iii) the line connecting P2 and Px = {qf, 0). 

If + ^ < 1, it is described by two lines: (i) the line connecting Py and P3 = {qf — Ai, — A2) and (ii) the 
line connecting P3 and Px- 

Proof: Please refer to |[T6l . ■ 

In Fig. 1 and Fig. 2, we illustrate the stability region of CARA with imperfect CSI along with the case with 
perfect CSI and the original ALOHA with no CSI. It is evident that the stabiUty region of CARA with perfect CSI 
always includes that of original ALOHA, which complies with the previous results obtained in |0- It is also obvious 
that the stability region of CARA with perfect CSI always includes that with imperfect CSI, and the difference 
between the two regions, therefore, can be understood as the loss due to the errors in channel estimation. However, 
the stability region of original ALOHA is not a proper subset of that of CARA with imperfect CSI as shown in the 
figures. The inefficiency of CARA is due to the fact that each node transmits only when the channel is estimated to 
be good and those time slots when the channel is estimated to be bad but when it is indeed good is not exploited. 



8 





Fig. 3. Stability region of CARA: the case of convex region (parameter setting:7rp = 0.8, tt^ = 0.7, qi\{G} = 1, ?2|{g} = 0.9, 

gi|{G,S} = 92|{S,G} = 0.5, 9l|{G,G} = g2|{G,G} = 0.4, £^ = 0.1,Vi, j) 
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IV. Analysis Using Stochastic Dominance 

In this section, we provide details on the derivation of our main result presented in the previous section. Note 
that nodes are interfering with each other only when they are transmitting and, under the considered protocol, node 
i transmits with probability pi if the estimated channel state is good and its queue is non-empty. Consequently, the 
service process of a queue at a node depends on the status of the queue at the other node and, thus, the content 
size of the queues form a two-dimensional discrete-time Markov chain with infinite state space, which makes the 
problem challenging even for the two-node case. Having to deal with a system with such interacting queues, a 
tool that can be used is the stochastic dominance technique which was introduced in |[T9l and further used in llT6l . 
|[2]l - |[23l . The essence of the stochastic dominance technique is to decouple the interaction between queues via the 
construction of a hypothetical system; this hypothetical system operates as follows: i) the packet arrivals at each 
node occur at exactly the same instants as in the original system, ii) the coin tosses that determine the transmission 
attempts of the nodes have exactly the same outcomes in both systems, iii) however, one of the nodes in the 
system continues to transmit dummy packets even when its packet queue is empty. Sending dummy packets is only 
aimed to cause constant interference to the other node and does not contribute to throughput if the transmission 
is successful. It is obvious that sample-pathwise the queue sizes in this dominant system will never be smaller 
than their counterparts in the original system, provided the queues start with identical initial conditions. Thus, the 
stability condition obtained for the dominant system is a sufficient condition for the stability of the original system. 
It turns out, however, that it is indeed sufficient and necessary, which will be discussed in detail later in this section. 

Construct a hypothetical system which is identical to the original system except that node 2 transmits dummy 
packets when it decides to transmit but its packet queue is empty. Define 1 j (n) to be an indicator function whose 
value is one if the transmission by node i is successful, which also necessarily requires that the corresponding 
node transmits at that time slot. Otherwise, li{n) = 0. By conditioning on the underlying actual channel states, the 
average service rate of the queue at node 1 can be expressed as 

= ^Pr[li|C = C']Pi[C = C] 
c 

= Pr[li|C = {G,G}]7rf7r2^ + Pr[li|C = {G, i?}]7rf vrf (6) 

where we used the fact that the transmission success probability is zero when a node's own channel state is bad 
and the time index n is suppressed in the steady-state. Note that each node i transmits with probability pi only 
when the estimated channel state is good. By further conditioning on the estimated channel states, it is expressed 
as 

Pr[li|C = C'] = Y^ Pr[li|C = C',C = c']Fr[C = C'\C = C] 
c' 

where 

Pr[li|C = {G, G}] = [(Zi|{G}Pi(l - V2) + gi|{G,G}PiP2] ef ef + (7) 
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and 

Pr[li|C = {G,5}] = (7i|{G}Piefef + [gi|{G}Pi(l " P2) + (/i|{b,g}PiP2] (8) 

By substituting ^ and © into ^ and, after some manipulations, the average service rate of the queue at node 1 
is derived as 

/"I = Trf efpi {qi\{G} - *!P2) 

where ^| is defined in Section Jill By Loynes' Theorem, the queue at node 1 is stable if Ai < ^1. Note that the 
queue size at node 1 in this dominant system can be modeled as a discrete-time M/M/l system with arrival rate 
Ai and the service rate given above, which does not depend on the status of the queue at node 2. For stable input 
rate Ai that is less than ^i, the queue at node 1 empties out with probability 

Pr[Q,=0] = l-^ = l 



Observe that the service process of the queue at node 2 depends on the status of the queue at node 1 since node 
1 is able to transmit only when its queue is non-empty and, thereby, interfering with node 2. By conditioning on 
the emptiness of the queue at node 1, the average service rate of the queue at node 2 can be expressed as 

= Pr[l2|Qi / 0]Pr[Qi / 0] + Pr[l2|Qi = 0]Pr[Qi = 0] (9) 

The service rate of the queue at node 2 when the queue at node 1 is non-empty can be obtained by following the 
same procedure used for deriving the service rate of the queue at node 1 and is given by 

Pr[li|Q2 / 0] = vrf efpi {q,\{G} " ^Pi) (10) 

The service rate of the queue at node 2 when the queue at node 1 is empty can be obtained quite simply as it does 
not depend on the action made by node 1 and is given by 

Pr[li|g2 = 0] =7rfe-fpigi||G} (H) 

Substituting dTOl) and (fTTl) into ^ yields the average service rate of the queue at node 2 which can be summarized 
to 

r-r ( ^2^1 \ 

and by Loyne's Theorem, the queue at node 2 is stable if A2 < /i2- Consequently, the stable input rate pairs (Ai, A2) 
are those less than (^1,^2) elementwise, and it gives the description of 7^2 in Lemma |3TT] By reversing the roles 
of the two nodes, we construct a parallel dominant system in which node 1 transmits dummy packets and, by 
following the same procedure used for the first dominant system, the stability region for this parallel dominant 
system is obtained as described in Lemma 13.11 which is denoted by TZi . 

Importantly, the stability condition obtained using the dominant system technique is not merely a sufficient 
condition for the stability of the original system but is sufficient and necessary. Consider, for example, the subregion 
IZ2 obtained for the dominant system in which node 2 transmits dummy packets. The sufficient part is trivial which 
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follows from the construction of the dominant system such that the dominant system stochastically dominates the 
original system in number of packets in the queues at all times. The necessary part can be proved as follows: if 
for some A2, the queue at node 2 is unstable in the hypothetical system, then Q2{n) approaches infinity almost 
surely. Note that as long as the queue does not empty, the behavior of the hypothetical system and the original 
system are identical, provided they start from the same initial conditions, since dummy packets will never have to 
be used. A sample -path that goes to infinity without visiting the empty state, which is a feasible one for a queue 
that is unstable, will be identical for both the hypothetical and the original systems. Therefore, the instability of 
the hypothetical system impUes the instability of the original system. 

We now obtain the maximum achievable stability region, which is the closure of the stability region 6(e,p) over 
the transmission probability vector p. An equivalent way of taking the closure operation is to optimize the boundary 
of the stabiHty region (S(e,p) over p. For the subregion Tl2, for example, we set up the following optimization 
problem, in which //2 is maximized over p at given Ai while satisfying the stability of the queue at node 1. 

max ,2 = .§^p2 - ^^O^J^^^^^^) (12) 

subject to Ai < vrf efpi (gi|{G} - ^1^2) (13) 
pe[0,lp (14) 

Note that 112 depends only on p2 but not on pi . Differentiating p2 with respect to p2 gives 

9^2 G-G ( ^2(ll\{G}^l \ 

By differentiating once again, we have 

_ 2^f^|7rQ^gi|{G}Ai 

dpi '^?^{qi\{G}-nP2f 

Observe that 

«1|{G} - ^lP2 > qi\{G} - ^1 

= (1 - TT^ef - )9i|{G} + TT^ef «i|{G,G} + Trf ef 

> (1 - - )gi|{G} + '^2 (-2 (ll\{G,G] + T^2 ^2 1l\{G,B} 

G —G B B 

= 7I"2 ^2 ?1|{G,G} + 7r2 62 qi\{G,B} 

> (15) 

Therefore, the second derivative is strictly negative and, consequently, p2 is a concave function of p2. Equating the 
first derivative to zero gives the maximizing P2 as 



*1 V V ^1 ^«2|{G} / 
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and the corresponding maximum function value is obtained by substituting ( [T6l ) into (fT2l ) as 



_ vr2^et / / ^igi|{G}Ai \ / / ^ig2|{G}Ai \ 

/^2,curve " 9l|{G} " \ n^G^ , ^ , 92|{G} " -l / G^„ , 

Suppose that the maximum occurs at a strictly interior point of the feasible region, i.e., P2 G (0, 1), which 
corresponds to the condition 

T^1^Q2\{G}{Qi\{G}-'^1? ^ . ^ T^l^qi\{G}Q2\{G} 

that is obtained by rearranging ([16} and substituting the extreme values of i.e., and 1. On the other hand, the 
constraint ([T3] ) should also be satisfied for the derived which gives 

A, < <!LflM£l (18) 

Q2\{G} 

Consequently, /ig ^urvc valid only for the intersection of (ITtT i and ([TST l. By comparing the endpoints, the intersection 
is specified to be the same with ([TT] ) if ^ 92\{G}^ which is impossible. This is because, from the relation described 
in (fTSl ). it can be deduced that "^f < qi\sr\- If — I — > 1, the intersection becomes 

* ^'ll^/ <?1|{G} <?2|{G} — 

7rfefg2|{G}('7i|{G}-^f)^ , Trfef^-lQiKG} 

re < Ai < 

^29l|{G} 92|{G} 

Otherwise, it is an empty set. 

Next suppose that is either or 1, which is the case when Ai is outside of the range ([TTl i. If Ai is on the 
right-hand side of the range, becomes non-positive and, thus, ^2 is a non-increasing function of p2. Therefore, 
P2 = 0, which gives ^2 = 0- Whereas, if Ai is on the left-hand side of the range, 1x2 is a non-decreasing function 
of p2 and, thus, P2 = 1. The corresponding maximum function value is obtained as 

* G-G ( ^2^1 \ 
/^2,line = -^2^2 92|{G} 

V ^1 ei (9i|{G} - ^!)/ 

By substituting P2 = ^ into (fT3l) . we have 

Ai < et((7i|{G} - ^1) (19) 

Given the fact that Ai lies on the left-hand side of the range ([TT] ) in addition to the above constraint, it is shown that 
if — ^ 1 2_ > 1 M* is valid for Ai on the entire range of the left-hand side of (iTTb. Otherwise, it is valid 

gi|{G} 12\{G) — ' ^^2,11110 J- 'I ' ' 

for the range specified by ([19). Following the similar procedure, we can optimize the boundary of the subregion 
IZi , which completes the proof of Theorem 13.11 

V. Comparisons with the Centralized Schedulers 

In this section, we compare the stability region of CARA to that achieved by scheduling policies that make 
centralized decision based on the CSI feedback. Note that a scheduler allocates each time slot to one of the nodes 
such that the scheduled node can transmit in an interference-free environment during the allocated slot. In [1311 . it 
was shown that queue length information can be utilized to improve the scheduling performance. Specifically, the 
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discovered policy that serves the longest-connected-queue (LCQ) among those in 'Good' channel state stabilizes 
the system whenever the input rate vector is inside the stability region. Here, the stability region is defined as the 
set of arrival rate vectors that can be stably supported by considering all possible stationary scheduling policies. 
This is why the LCQ policy is called a throughput-optimal policy. The following theorem is derived again for the 
case with channel estimation errors, which was originally derived in ||3TI for the case with perfect CSI. 

Theorem 5.1: The necessary and sufficient stability condition by considering all possible stationary scheduling 
policies in the presence of channel estimation errors is 

j:^<i- n vAA'G{i,...,iv} (20) 

Furthermore, the LCQ policy stabilizes the system as long as it is stabilizable. 

Proof: Assume that the system is operating under certain stationary policy and it is stable. Denote by lj(n) the 
indicator function that is equal to 1 if the transmission by node i is successful, which necessarily implies that node 
i is chosen to transmit at that time slot. Also denote by h (n) the indicator function that is equal to 1 if the actual 
channel state between node i and the receiver is Good and it is estimated correctly. The expectations are given by 
= qi\{c^{n)}^ ^nd E[Ii{n)\ = T^fef. The number of packets at queue i evolves with time according to the 
queueing dynamics in equation ([T) with ^i{n) = hi{n) ■ li{n), where hi{t) = 1 if node i is scheduled at time slot 
n and Ii{n) = 1. Thus, the departure rate from queue i is written as 

E[h^{t)li{t)] = qi\{G}E[hi{t)\ 

Note that if the system is stable, the rate of what comes in must be equal to the rate of what goes out. In other 
words, for any subset of network nodes M' , the following equality must hold. 

E ^ = E ^[^'(*)] (21) 

Consider now the partition of the probability space into the events 

Jli = {Ii{t) =0,i G AA'} 

1^2 = { /i (t) = 0, i e AA' }^ n { Qi(t - 1) = 0, i G AA' } 
1^3 = {/i(t) =0,i GAA'fn{Qi(i-l) = 0,i G AA'f 
where U,^ is the complementary set of Q. Notice that E [X^igAT' = for / = 1, 2 and, therefore, we have 



Y,E[hi{t)] = E 



E 



< 1 - Pr[J^i] - Pr[Jl2] 

Owing to the assumption on the independence of the channel processes between different nodes, we have 



Pr[J7i] = n (1 
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and Pr[il2] > 0. Therefore, we have 

Y,EMt)]<l - [](l-vrf6f) (22) 

ieAA' ieAf' 

and equations (|2TI) and (|22]) implies equation (l20l ) in the theorem. 

For the sufficiency, we use the Lyapunov drift argument similar to that used in IIBTI . Denote by L{Q{t)) a 
Lyapunov function of the queue length process and define the conditional Lyapunov drift as A{L{Q{t))) = 
E[L{Q{t + 1)) — L{Q(t))\Q{t)]. If there exist some e > and a finite number b such that the conditional Lyapunov 
drift satisfies A{L{Q{t))) < —e for L{Q{t)) > b, the queues in the system are stable. This is an application of 
Foster's criterion for ergodicity of a Markov chain [27 1. Let us consider L{Q) = Xli^i liunyQi - Then, for the 



considered Lyapunov function, it can be shown that the conditional Lyapunov drift satisfies 



N 



A(L(g(t))) < 1 + J] qJa} + 2A,Q,,(t)) - 2E 



1=1 



N 



^Qi{t)h,{t + l)\Qit) 



.i=l 



(23) 



The details on the derivation of the above inequality can be found in |[3T| . 

Define now a permutation for z = 0, . . . , such that eo = and (i) > Qei_i(i)> for i = 2, . 
Consider also a partition of the probability space into the events defined by 

= {IeAt + 1) = + 1) = 0, N>j>i} 

for i = I, . . . , N and $o = {-^(^ + 1) = 0}. Notice that from the definition of the LCQ policy, in the event <J>j, 
queue Cj is served if it is not empty. Therefore, the last term in the right-hand side of ( [231 ) becomes 



E 



N 



Y,Q^{t)hi{t + l)\Q{t) 



1=1 



By substituting the above into 



N N 
2=1 j=i+l 

and after some manipulation, we obtain 

:{i,...,7V} 

It is not difficult to observe that if the condition ([20] ) is met, then for sufficiently large Qej^{t), the right-hand side 
of the above conditional Lyapunov drift becomes negative. This completes the proof. ■ 



N 



A{L{Qm <1 + Y1 E[AUt)] + 2Q,At) max 
'1=1 



f) 



For the two-node case, the stability condition ([20] | in the above theorem becomes 



Ai 



+ 



A, 



< TTi ei + £2 (1 - ^1 ei ) 



G-Gi 



1l\{G} <i2\{G} 

and Aj < '^f^(li\{G} for ^ ^ {li^}. In Fig. 2, we compare the stability region of CARA to that achieved by 
LCQ policy, and it can be observed that the former is not necessarily a subset of the latter. Indeed, the relationship 
between them depends on parameters such as the channel estimation error and the MPR probabilities. Note that 
the stability region of CARA becomes a subset of that of the LCQ policy when the stability region of CARA is 
non-convex, or convex but if the vertex P3 given in Q is strictly contained in the stability region achieved by LCQ 
policy. The condition for i-3 to be inside the stability region of the LCQ policy is given by 



—C 

^2 ^2 ?1|{G} 



+ 



T^l^<i2\{G} 



> 1 
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— □ — LCQ policy (parameter setting 1 ) 
CARA (parameter setting 1 ) 

- B - LCQ policy (parameter setting 2) 

- ^ - CARA (parameter setting 2) 



0.8 - 
0.7 - 
0.6 - 




Fig. 4. Comparison with the LCQ policy (common setting: Trf = 0.8, 71-2' = 0.7, qi\{G} = 92|{g} = 0.9, parameter setting 1: qi\{G,B} = 
q2\{B,G} = 0.7, gi|{G,G} = g2|{G,G} = 0.6, €^ = 0.1,Vi, j, parameter setting 2: gi|{G,s} = 92|{s,g} = 0.4, gi|{G,G} = 92|{g,g} = 0.3, 
ef =0.3,Vi,j) 

Otherwise, if the stability region of CARA is convex and the above inequality does not hold, the stability region 
of CARA is not a proper subset of that of the LCQ policy, i.e., there exists a region that can be achieved only by 
CARA. 

VI. Concluding Remarks 

In this work, we studied the stability property of CARA in the presence of channel estimation errors and showed 
that its stabiUty region may not strictly contain that of original ALOHA. To guarantee the superiority of CARA 
even with imperfect CSI, we need to modify the protocol itself such that each node transmits with some positive 
probability, although it believes that the channel is in the bad state. Such modification was not considered here. We 

also compared the stability region of CARA to that achieved by the throughput-optimal LCQ policy and showed 
that the former is not necessarily a subset of the latter especially as the MPR capability improves. The stability 
analysis of CARA had to resolve the complex interaction between nodes and, hence, extending the results to the 
general case with an arbitrary number of nodes, although highly desirable, presents serious difficulties and was not 
considered here. 
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